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Abstract 

We simulate two dimensional QED with two degenerate Wilson fermions and plaquette 
gauge action. As a consequence of the Mermin- Wagner theorem, in the continuum limit 
chiral symmetry is realized a la Wigner. This property affects also the size of the cutoff 
effects. That can be understood in view of the fact that the leading lattice artifacts are 
described, in the continuum Symanzik effective theory, by chirality breaking terms. In 
particular, vacuum expectation values of non-chirality-breaking operators are expected 
to be 0(a) improved in the chiral limit. We provide a numerical confirmation of this 
expectation by performing a scaling test. 
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1. Cutoff effects in lattice (gauge) theories can be described using an effective 
continuum action, as proposed by Symanzik in refs. [1,2]. In this approach the 
leading lattice artifacts (e.g. in the spectrum of the theory) can be removed by 
including a set of irrelevant operators in the action and by properly tuning their 
coefficients. For the case of the Wilson lattice regularization of QCD [3], the 
relevant coefficient can be tuned by requiring the restoration of chiral symmetry 
up to O(a^). This interplay between chiral symmetry and cutoff effects has been 
addressed in detail in ref. [4]. 

Further insights on this connection have been recently derived in ref. [5] by 
considering also so called spurionic lattice symmetries to classify the operators, 
which can appear in the Symanzik effective theory. Without reproducing the 
whole argument, here we will simply summarize the results relevant as premises 
for this work. Let's consider the Wilson fermionic action for two degenerate ffavors 

Sf[U,iP,^] = a'^''^^ip(x){D + mo)ip(x) and (1) 

^ = H7m(V; + - arV;v4 , (2) 

where V* and V,, are the covariant backward and forward lattice derivative re- 
spectively, U denotes the gauge field and r is the Wilson parameter, which we set 
to 1. The vacuum expectation value of a multiplicatively renormalizable operator 
O can be expanded as 

mr,m, = [C^ + am,e] mZ + aJ^i^^rvS^O^Z + ' (3) 

I 

where uiq is the bare subtracted fermion mass defined as rriq = rriQ — nic, such 
that the physical fermion mass vanishes for mo = rric. We refer to [5] for any 
unexplained notation in eq. Q. The operators Oi appearing on the rhs result 
from the insertion of the 0(a) terms in the action and from the 0(a) terms 
associated with the operator O itself. They can be classified according to their 
parity Pj^^ under the transformation 

i?5 : V ^ = 75^^ , ifj -^ip' = -tjj-f^ , (4) 

which is a non-anomalous element of the chiral group and produces a spurionic 
symmetry of the Wilson action when combined with the replacements r — *• — r, 
TTiq — i> —rriq. The authors of ref. [5] have shown that 

Pg + PS:+ni = I mod (2), (5) 

which, loosely speaking, implies that the 0(a) terms in the chiral limit (where the 
sum reduces to the ni = contributions) have opposite i?5-parity compared to 



1 



the leading term. It is indeed interesting to consider the hmit rriq — > in eq. (jHl). 
Two different scenarios are possible 

• Spontaneous breaking of chiral symmetry does not occur, as for QCD in 
small volume, therefore the theory is analytical at rriq = and we can 
directly set the fermion mass to zero in eq. (jH)). From that we infer that if 
O is even under then the operators Oi are odd according to eq. (0) and 
their vacuum expectation values vanish in the continuum (because of chiral 
symmetry). We conclude that in this case (0)|r,mq is free from 0(a) effects 
in the chiral limit. Conversely, if P^^ = 1 mod (2) the continuum limit of its 
vacuum expectation value (vanishing for symmetry reasons) is approached 
with a rate proportional to a. 

• Chiral symmetry is realized a la Goldstone. In this case, due to the non- 
analyticity at rriq = 0, the chiral point can only be approached through a 
limiting procedure. Still, "automatically" 0(a) improved correlation func- 
tions can be obtained using Wilson- or mass-averages or, more practically, 
by employing twisted mass fermions at maximal twist [5-7]. 

All these considerations apply to any fermionic theory regularized a la Wilson. 
In particular we want to numerically test the first scenario described above by 
considering the Schwinger model [8] with two dynamical flavors, such that the 
transformation is well defined. More importantly, in two dimensions continuous 
chiral symmetry cannot be spontaneously broken due to the Mermin- Wagner the- 
orem [9] . Unfortunately, mainly for numerical reasons, we will not be able to work 
with massless fermions. Therefore in addition to the O(a^) cutoff effects expected 
in the chiral limit we might observe O(amq) effects on our quantities. 

2. We simulated two dimensional compact QED on a torus with periodic boundary 
conditions in time and space. Since the gauge coupling g is of mass dimension 
one, the model is super-renormalizable. For the lattice theory this implies that the 
continuum limit in a finite physical volume can be taken at fixed g-L (see ref. [10]). 
For later usage we introduce the dimensionless coupling f3 = (ag)^"^. Clearly, when 
taking the continuum limit a suitable fermion mass m has to be kept fixed as well. 
We decided to define m through the PCAC relation [4, 11, 12] (see also below for 
details) and fixed the product m ■ L to a constant value. Notice that due to the 
super-renormalizability of the model we do not need to compute renormalization 
factors Z and can just use the bare PCAC mass. Indeed, in perturbation theory 
Z can be written as Z = 1 + Z^^'^a^g'^ + . . . , therefore, at fixed g, loop corrections 
only change the O(a^) ambiguities. Similarly, if we wanted to fully 0(a) improve 
the theory (and remove 0(a) effects also from vacuum expectation values of -R5 
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odd operators) by adding the Sheikholeslami-Wohlert term [13] to the action, its 
coefficient could be set to 1 to all orders in the perturbative expansion. The same 
is true for the 0(a) counterterms of the operators. In other words the 0(a) cutoff 
effects, if any, are tree level cutoff effects. 

For the simulations we used the Hybrid Monte Carlo algorithm [14] with a 
leapfrog integration scheme. Observables are constructed from the correlation 
functions (we set a = 1 and write x = (xo,x)) 

Cxy(xo) = E ^(0' ^) n ^)rxrXO, y) ^{xo, z)rYrXxo, z) , (6) 

x,y,z i=x 

with X and Y = A or P and 

Ta = 7o75 , Tp = 75 , (7) 

while the Pauli matrices r" act on flavor indices. In eq. © the product of space- 
like gauge links is needed to deflne gauge invariant wall-to-wall correlators. The 
additional numerical effort required to construct such correlation functions is quite 
moderate in two dimensions. 

The PCAC mass m is computed through the ratio 

= m , (8) 

2Cpp(xo) 

derived from the axial Ward identity. Our scaling quantities are obtained from 
the correlator Caa{xo), which for xq around T/2 is expected to be dominated by 
the lowest zero momentum state vr in the pseudoscalar sector. In this case, the 
correlator is described by 

Caa(xo) = <I>1 cosh(m^(r/2 - xq)) , for xq ~ T/2 , (9) 

where the cosh function is due to the periodicity in time and the matrix element 
$^ is, up to the normalization, the analogon of the pion decay constant in 
QCD. We will see in the next section that the formula in eq. Q reproduces the 
data fairly well, which is plausible as we have m,r-^ — 5 and T = 2L. Since the 
correlator Caa{xo) is clearly even under we expect the dimensionless quantities 
LitIt^ and L^t^ to approach their continuum limit values with a rate proportional 
to a^ up to corrections of O(amq). 

3. The simulation parameters are collected in table Q together with the results. 
We also give some details concerning the algorithm. The hopping parameter 
K = (2mo + 4)~^ is tuned in order the keep the fermion mass constant within 
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L/a 







rriL 




TV 




TT J 


accept. 


16 


2 


0.2680 


1.01(1) 


4.80(6) 


0.0400(7) 


50 


5000 


96% 


20 


3.125 


0.2603 


1.00(1) 


4.7(1) 


0.035(1) 


50 


5000 


94% 


24 


4.5 


0.2564 


1.008(7) 


4.7(1) 


0.0321(6) 


50 


4000 


93% 


32 


8 


0.2530 


0.995(6) 


4.7(1) 


0.0276(9) 


60 


2500 


93% 


40 


12.5 


0.25153 


1.004(8) 


4.68(8) 


0.0269(9) 


70 


1500 


92% 



Table 1: Simulation parameters and results. The length of the single trajectory is 
always 1, discretized in ngtep intervals. 



1% accuracy. To extract the pseudoscalar mass we define a local effective mass, 
which assumes the correlator CAA(a;o) to be dominated by a single state, and 
we average it over a plateau region. Similarly, we compute by averaging the 
ratio [Caa{xo)/ cosh{mT,{T /2 — xq))]^^"^ over the same region. The results for Lm^r 
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Figure 1: Scaling plot for the pseudoscalar mass m^r. 

are plotted in figure ^ against a/L. It is clear from the plot that within the 2% 
errors we do not see any cutoff effect. The symbols at a/L = correspond to the 
predictions obtained for our value of the fermion mass from different approximate 
analytical solutions valid in the limit of small mass and large coupling g [15]. We 
regard the observed consistency as a check of our setup. In addition, for the same 
quantity and for a similar choice of parameters, results consistent with lattice 
artifacts linear in have been recently reported also in ref. [16]. 

The discussion of is a bit more delicate since we see cutoff effects in this 
quantity. As it is shown in figure El those are clearly consistent with being linear in 
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Figure 2: Scaling plot for the matrix element $7r- 

o} only. Nevertheless, in order to estimate the size of the O(amq) effects, we tried 
to fit the data also to a polynomial with terms linear and quadratic in a. The fit is 
acceptable in terms of and the continuum limit we obtain is in agreement with 
the one in figure |2l but it has a five times larger error. The coefficients of the linear 
and quadratic terms have large errors as well. They are both consistent with zero 
but strongly anticorrelated. We conclude that the sensitivity of our data to the 
O(amq) effects is very small. Adding a smaller lattice resolution would probably 
help to disentangle them from the O(a^). 

4. The numerical study presented here confirms the expectation that in the 
absence of spontaneous chiral symmetry breaking cutoff effects are of O(a^) also 
when Wilson fermions are used (at least if the chiral limit is considered). The 
situation is very different from QCD in four dimensions (and large volume), where 
for Wilson fermions the 0(a) effects are rather large [17] and have to be removed 
by following the Symanzik improvement programme [4]. 

As a consequence, testing fermionic actions by scaling studies in the Schwinger 
model provides, in our opinion, very little information about the cutoff effects for 
the same regularizations in the phenomenologically more relevant case of QCD. 

On the other hand, to improve our confidence in the argument presented 
here, it would be interesting to extend the study by considering different values of 
the fermion mass in order to assess more precisely the size of the residual O(amq) 
effects. As far as we can tell now, those appear to be fairly small. In addition, the 
mass of the scalar particle r] could be included among the observables. Contrary to 
the quantities discussed here, this mass does not vanish in the chiral limit [12]. To 
this end, the numerical techniques introduced in ref. [18] could provide an efficient 
way to evaluate the contributions coming from disconnected quark diagrams. 
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